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INTRODUCCION 


Introduction - Biblioarafia 


Esta publicacion pretende brindar en una forma rapida y sencilla los conceptos basicos de 
la dinamica de estructuras aplicada a las construcciones civiles, enfocada desde un punto 
de vista numerico. Su organizacion en capftulos se basa en la secuencia de unidades tema- 
ticas del programa de la asignatura homonima que se dicta en la carrera de Ingenierfa Civil 
de la Facultad Regional Parana de la Universidad Tecnologica Nacional. 

Todos los temas se desarrollan como sintesis de libros y publicaciones de otros autores 
(ver Bibliograffa) solo que organizados segun la conveniencia que mejor se ajusta al curso 
impartido. 

El Capitulo 1 hace una muy breve resena a conceptos sobre sismologfa y su relacion con la 
ingenierfa estructural, caracterizacion, registros, etc. 

El Capitulo 2 brinda una aproximacion a los conceptos basicos de la dinamica estructural, 
la importancia de la masa, la relacion entre velocidad de carga y de reaccion de una estruc- 
tura. Avanza sobre los modelos estructurales dinamicos, grados de libertad y la discretiza- 
cion espacial y temporal. 

El Capitulo 3 esta dedicado a la caracterizacion y calculo de osciladores de un grado de 
libertad dinamico. Son tratadas las vibraciones libres y forzadas, estas ultimas con cargas 
armonicas y arbitrarias. Se presenta la resolucion numerica mediante integracion directa. 

El Capitulo 4 esta dedicado a las estructuras con multiples grados de libertad dinamicos. 
Se estudian las vibraciones libres y el calculo de modos y frecuencias propias, tambien se 
presenta la resolucion de los sistemas dinamicos mediante diversos metodos como son: 
descomposicion y superposicion modal, integracion directa y respuesta maxima mediante 
espectros de respuesta. 

En el Capitulo 5 se realiza una aproximacion al tratamiento dinamico de los efectos del 
viento sobre las construcciones teniendo como base el reglamento CIRSOC 102. 

Por ultimo en el Capitulo 6 son tratados los efectos sfsmicos y el analisis estructural, desde 
el punto de la dinamica estructural, enfatizando en la definicion de la accion sfsmica y la 
comparacion de las formulas vistas en capftulos anteriores con los reglamentos vigentes. 


Un agmdecimiento especial a Marlene Jaurena, que con paciencia y dedi- 
cacion realizo la transcripcion de los manuscritos de este trabajo. 
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Capitulo 1 - Nociones Basicas de Sismologia 


CAPITULO 1 - NOCIONES BASICAS DE SISMOLOGIA 
Causas que generan los terremotos o sismos 

Los terremotos pueden definirse como movimientos de la corteza terrestre, con amplitudes 
y frecuencias dependientes del tiempo. 

Las causas que los generan son variadas: 

Terremotos de colapso : son los originados en cavidades subterraneas por el colapso de las 
mismas, son de baja intensidad. 

Terremotos de origen volcanico: la explosion de gases durante las erupciones volcanicas 
puede producir terremotos que, en general, tienen una intensidad pequena y afectan a su¬ 
perficies limitadas. 

Terremotos tectonicos: estan causados por la rotura brusca de las capas rocosas a lo largo 
de superficies de fractura (falias), son los mas fuertes y mas frecuentes. 

Terremotos causados por explosiones: las explosiones producidas por el hombre son capa- 
ces de generar vibraciones del terreno, con una intensidad tal que pueda causar movimien¬ 
tos en las estructuras. 

En general, el movimiento de la corteza se produce por un choque o movimiento brusco 
ocurrido a una cierta profundidad bajo la superficie terrestre en un punto teorico denomi- 
nado foco o hipocentro, a su proyeccion sobre la superficie terrestre se le denomina epi- 
centro. 



Fig. 1.1 - Definiciones geometricas de un sismo 
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Los sismos desde el punto de vista de la ingenierfa y su caracteri- 
zacion 

Los terremotos mas importantes son los tectonicos, pues son los que traen consecuencias 
mas desastrosas en las estracturas que afectan, debido a esto, son los que se tienen en 
cuenta para la elaboracion de normas para la contraccion de estructuras sismoresistentes. 


La intensidad sismica es una medida de los efectos de los terremotos en el entorno y en 
particular sobre las estructuras. 

Existen diferentes escalas de intensidades que describen, para cada valor que esta tome, 
los efectos que produce el terremoto. Una de las mas difundidas es la escala de Mercalli 
Modificada. 

Algunos de los efectos sobre las estructuras en orden creciente de intensidad son: 

1. fisuracion de las estructuras de madera 

2. agrietamiento de las estructuras debiles de mampostena 

3. agrietamiento de las estructuras ordinarias de mampostena 

4. colapso parcial de estructuras ordinarias de mampostena; dano en estructuras bien 
ejecutadas de mampostena no disenadas para resistir fuerzas sismicas 

5. colapso de estructuras ordinarias de mampostena; las estructuras con diseno anti- 
sismico son seriamente danadas; danos en cimientos; grietas en el terreno 

la mayorfa de las estructuras son destruidas junto con sus cimientos, danos importantes en 
presas y diques, grandes deslizamientos del terreno destruccion casi total, grandes masas 
de rocas desplazadas, etc. 

Un sismo se caracteriza por su intensidad (parametro subjetivo) y por su magnitud (pa¬ 
rametro objetivo). 

La escala objetiva mas popular es la de Ritcher, en la que la magnitud M mide la energia 
del terremoto en el foco y es el logaritmo decimal de la amplitud del movimiento sismico, 
medido en micrones a 100[km] del epicentro, por un sismografo Wood-Anderson estan- 
dar. La magnitud M esta relacionada con la energia del terremoto, en ergios, por la expre- 
sion: 

log E = 1 1,8 + 1,5M 

Se han establecido varias relaciones empiricas entre la intensidad I MM y la magnitud M, 
enumeramos algunas a continuacion: 

Esteva y Rosenblueth: 

I MM = 8,16 +1,45 M - 2,46 log R 
R: distancia focal en [km] 

Tambien se ha relacionado la magnitud M con los valores maximos de las caracteristicas 
cinematicas del movimiento, estas relaciones se han establecido estadisticamente. 
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Donovan: 


_ 1080 e°’ 5M 
a '" ~ lR + 25 f 25 

a m : aceleracion maxima del terreno en [cm/s 2 ] 
R: distancia focal en [km] 

Esteva y Villaverde: 

C ^7 0,8 M 

<*m _ 5,7 g 
8 {R + 40) 2 

32 e M 

V " ~ (JJ + 25) 1 ’ 7 

a m : aceleracion maxima en [cm/s 2 ] 
v m : velocidad maxima en [cm/s] 

R: distancia focal en [km] 

Esteva: 

a m =1230e°’ 8M (R + 25)' 2 

v m = 15 e M (r + 0,17 e 0,59M )~ 17 

a m : aceleracion maxima en [cm/s'] 
v m : velocidad maxima en [cm/s] 

R: distancia focal en [km] 


Registro de ondas si'smicas. Parametros utilizados y mapas de 
riesgo si'smico 

Los terremotos son fenomenos debidos a la brusca liberacion de la energfa de deformacion 
acumulada durante largos periodos de tiempo en la zona superficial de la tierra. Los sis- 
mos producen ondas de varios tipos, que se propagan desde su foco en todas las direccio- 
nes a traves de la tierra. Estas ondas son registradas mediante aparatos denominados sis- 
mografos, disenados para medir la aceleracion, la velocidad o el desplazamiento del mo- 
vimiento sfsmico. Estos parametros son relativos, ya que los valores obtenidos estan afec- 
tados por las caracterfsticas del instrumento registrador y por las condiciones de ruido am- 
biental en el lugar de registro. 

Los mapas de riesgo sfsmico representan una sfntesis de todos los datos sismologicos y 
geologicos de un pats. Estos mapas se utilizan para determinar el nivel de proteccion que 
se debe alcanzar en las estructuras en cada zona de riesgo. 

Diversos aspectos brindan la subdivision en zonas, pero los fundamentales son: 
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Estudios geologicos y geotecnicos : proporcionan datos de composicion y caracteristicas 
dinamicas de las rocas y capas de suelo que componen la corteza terrestre. 

Estudios sismologicos : sintetizan los parametros que caracterizan la sismicidad de la zona: 

1. ubicacion de falias 

2. registro de los terremotos que ocurren en la zona 

3. mapas de epicentros 

4. datos histdricos 

5. periodos de retorno (intervalo medio de tiempo en que se espera ocurran dos sismos de 
igual o mayor intensidad) 

6. datos del mecanismo focal 

7. correlacion de la sismicidad de la zona analizada con la de la macrozona en la que se 
encuentra 

Estudios de Ingenieria y Sismologia : 

1. analisis del efecto que han producido sobre las estructuras y las personas los terremotos 
ocurridos en el pasado 

2. “prediccion ” estadistica de las caracteristicas mas probables de la accion sismica que 
se produzca en la zona 

Es importante destacar que la geologfa local de la zona puede modificar la propagacion de 
las ondas sfsmicas. Las ondas se reflejan y se refractan cuando en su recorrido aparece una 
discontinuidad, por ejemplo una variacion de las caracteristicas mecanicas del terreno, ello 
produce cambios en la velocidad. 

En general, el calculo y la cuantificacion de las acciones sfsmicas en la estructuras se rea- 
liza en funcion de protocolos, secuencias y definiciones de acciones dadas por normas y 
reglamentos. En los capftulos siguientes se ofrecen algunas aplicaciones de este tipo. 
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CAPITULO 2 - CONCEPTOS BASICOS DE DINAMICA ES- 

TRUCTURAL 

En un sentido amplio, un sistema dinamico es aquel cuyas variables experimentan varia- 
ciones en el tiempo y, si se conocen las influencias externas que actuan sobre el sistema, 
podra predecirse el comportamiento de este. 

SISTEMA variables con variaciones 

DINAMICO temporales 

conocidas estas acciones 
externas, permiten 
"predecir" el comportamiento 
de las variables temporales 



En nuestro curso, los sistemas a estudiar seran sistemas estructurales, las variaciones en el 
tiempo seran vibraciones producidas por cargas dinamicas. 


SISTEMAS 

ESTRUCTURALES 


cargas dinamicas 


permiten evaluar el 
comportamiento de 
la estructura frente 
a acciones dinamicas 




Definicion de la accion dinamica 

Una accion tiene caracter dinamico cuando su variacion con el tiempo es rapida y da ori- 
gen a fuerzas de inercia comparables en magnitud con las fuerzas estaticas. Algunas fuen- 
tes importantes de vibraciones estructurales son: 
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- sismos 

- viento 

- olas y corrientes de agua 

- explosiones e impactos 

- cargas moviles (vehiculos, personas, etc.) 


La definicion de estas cargas extemas puede distinguirse entre: determinista y no determi- 
nista, esta ultima denominada tambien estocastica o aleatoria. 

determinista : cuando su variacion temporal es perfectamente conocida 

no determinista : cuando alguno o todos sus pardmetros son definidos estadisticamente 

En nuestro curso trabajaremos con cargas definidas en forma DETERMINISTA. 


Respuesta dinamica —» cualquier magnitud que pueda caracterizar el efecto de una cargo 

dinamica sobre la estructura 


Una carga definida determimsticamente da origen a una respuesta, tambien determinista. 



Fig. 2.1- Definicion de la respuesta dinamica: para un punto considerado se calculan: 

deformaciones, aceleraciones, tensiones, etc. 


Acciones y fuerzas dinamicas 

Las acciones dinamicas definidas utilizando representaciones deterministas, son funciones 
del tiempo cuyo valor en cada instante ES CONOCIDO. 

Este tipo de representacion es apropiado para evaluar el comportamiento de una estructura 
A POSTERIORI del acontecimiento que dio lugar a dicha accion. Por ejemplo, evaluar el 
comportamiento de un edificio nuevo ante el terremoto ocurrido en Mexico en 1986 (del 
que se poseen registros). El diseno de una estructura NO PUEDE encararse en base a ac¬ 
ciones deterministas, pues nada nos asegura que la accion estudiada volvera a repetirse. 
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ACCION 


1 




este esquema temporal 
de carga debe ser 
perfectamente conocido 


Fig. 2.2 - Action y respuesta determinista 


ESTRUCTURA 


masa M 
rigidez K 


la comparacion se 
basa en: - velocidad de 
la accion 

ACCI0N " - peiiodo propio 

DINAMICA de la estructura 

da origen a fuerzas 

de "inercia" comparables j 

con las estaticas ' 

depende de 

K , M 


con las estaticas 


\ No considerada como 
"carga" sino como una 
propiedad intrfnseca de 
la estructura 

Fig. 2.3 - Action dindmica y propiedades de la estructura 


Importancia de la masa en el problema dinamico 

Aunque la carga varfe con el tiempo, la respuesta de una estructura varfa radicalmente se- 
gun la masa que vibra con ella. Ante una misma funcion de carga, una estructura SIN 
MASA y una CON MASA responden de la siguiente manera: 
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a) Estructura SIN MASA => SIN INERCIA 


x(t) 



masa: m = 0 



b) Estructura CON MASA => CON INERCIA 



masa: m 



la respuesta y deja vibraciones 


remanentes 


Fig. 2.4 - Importancia de la masa en la respuesta 
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Velocidad de reaccion de una estructura 

Ante una accion exterior, distintas estructuras reaccionaran de formas diferentes. Esta res- 
puesta esta mtimamente relacionada con las formas o modos de vibrar y sus correspon- 
dientes frecuencias o periodos propios. En el caso de un oscilador de 1 grado de libertad, 
este periodo propio se obtiene facilmente. No asf para estructuras de multiples GLD. 

Como veremos en los capitulos siguientes, los periodos y formas de vibrar dependen de las 
caracterfsticas geometricas y de materiales (rigidez) y de la inercia que la estructura opone 
al movimiento (masa). 



^ x 


m 



k 


/ 


/ 

w 

// 

t 

// 


Xo 


vibraciones 

libres 



k 


con AMORTIGUAMIENTO 
la amplitud decrece en cada ciclo 



practicamente constante 


F 


A 


tD 




En general si t D » T => no es necesario un analisis dinamico 
si to = T => PROBLEMA DINAMICO 


Fig. 2.5 - Velocidad de reaccion T vs. t D 
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m;k; 



a masa constante 
ki > k 2 

a rigidez constante 
m2 > mi 
t 


m, I m 7 I 

variables —> > —» T. ~ > —> , etc 

kj ‘ kj 

Fig. 2.6 - Velocidad de reaction; varios T, 


Modelos dinamicos caracteri'sticos 

Desde el punto de vista del calculo numerico, obtener la respuesta dinamica de una estruc- 
tura, es el resultado de "filtrar" la senal de excitacion a traves de la misma estructura y ob¬ 
tener las variaciones de las magnitudes de analisis (desplazamientos, velocidades, acelera- 
ciones, momentos, tensiones, etc.) respecto del tiempo. 

La obtencion de la respuesta requiere, previamente, la definition del movimiento del te- 
rreno (en caso sismico) tanto como de las caracterfsticas estructurales del mismo y de la 
estructura propiamente dicha. El analisis es practicado, no a la propia estructura sino a un 
modelo mecanico de la misma. La definition del modelo depende del tipo de estructura 
analizado y pretende brindar una serie de relaciones entre acciones y respuesta que descri- 
ban un modelo matematico del problema. 

Este modelo matematico puede ser resuelto mediante diversas tecnicas. En nuestro caso 
haremos hincapie en los metodos numericos de analisis. 

Segun la certeza con que fueron formulados los modelos y procedimientos o algoritmos de 
calculo durante el analisis, sera la precision de la respuesta obtenida. 
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original reflejada 


\ 

f\\J 

(\ f n 


1 

1 1 j 

V t 

onda sismica 


original 



Fig. 2.7 - “Filtmdo ” de una serial sismica 

Se brindan, a continuation, algunas definiciones tipicas del analisis estructural dinamico 
de una estructura: 


Grados de libertad (GL) 

Se definen como grados de libertad (GL) a los puntos de la estructura en los cuales se 
identifica algun desplazamiento y permiten brindar una deformada de la estructura. 


Grados de libertad dinamicos (GLD) 

Son los grados de libertad que tienen asociada masa y para los cuales puede conocerse las 
vibraciones o movimientos a lo largo del tiempo. 
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RESPUESTA (desplaz. en el piso 1) 


Fig. 2.8- Modelizacion de una estructura 


Metodos de modelizacion dinamica 

Pueden distinguirse modelos dinamicos exactos y modelos dinamicos discretos . 

En general, para la primera clase, solo pueden resolverse casos muy sencillos y con poca 
aplicacion practica, por lo que a lo largo del curso profundizaremos en modelos discretos. 

Para estos metodos modelos discretos, se debe tener en cuenta que la subdivision en domi- 
nios finitos es tanto espacial (discretizacion estmctural) como temporal (solucion para ins¬ 
tates de tiempo determinados). 
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n° infinito de 
puntos espaciales 

n° imfinito de 
puntos temporales 


DISCRETIZACION 

ESPACIAL 


DISCRETIZACION 

TEMPORAL 


Fig. 2.9 - Modelos dindmicos 

Discretizacion espacial de las estructuras 

Fundamentalmente, la diferencia con lo visto en otros cursos de analisis estructural (estati- 
co) radica en que en dinamica estructural, cuando hablamos de discretizar espacialmente, 
nos referimos a los GLD. 

Un modelo dinamico exacto (con infinitos GLD) acarrearfa mas inconvenientes en la reso- 
lucion matematica que beneficios en su precision. Ademas, en estructuras de edificios y en 
la mayorfa de las estructuras civiles, las masas se encuentran mas o menos concentradas en 
lugares conocidos. Es por esto que nuestro principal metodo de modelizacion dinamica 
sera el de las MASAS CONCENTRADAS . 

No obstante, existen otros, como ser: 

- metodo de los DESPLAZAMIENTOS GENERALIZADOS 

- metodo de los ELEMENTOS FINITOS 
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Metodo de las masas concentradas 



n° total de componentes 
de desplazamiento segun 
los cuales vibran las masas 
concentradas 


se distinguen: 

- MODELOS DE 1 GLD 

- MODELOS DE 
MULTIPLES GLD 


se calcula la "deformada" 
del modelo en cada instante 

V 

tensiones y deformaciones 
especificas pueden conocerse 
mediante los procedimientos 
del analisis estatico 


Modelos con 1 GLD: 

mm m 



(a) (b) v (c) 

Fig. 2.10 - Modelos con un solo grado de libertad. (a) modelo conservativo; (b) modelo 

con amortiguamiento; (c) modelo si'smico. 


m 



Fig. 2.11 - Estructuras modelizadas como un sistema de un solo grado de libertad. (a) 
portico; (b) el mismo portico con la masa concentrada al nivel de la viga; (c) modelo di- 

ndmico. 
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Modelos con multiples GLD: 



Fig. 2.12 -Modelos con varios grados de libertad. (a) modelo conservativo; (b) modelo 

con amortiguamiento; (c) modelo sismico. 



Fig. 2.13 - Estructura con dos grados de libertad: Portico de dos pisos y su modelo dind- 

mico. 



Fig. 2.14 - Estructura con masa distribuida (antena) y su modelo dindmico discreto con n 

grados de libertad. 
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Fig. 2.15 - Modelo dindmico de un portico de cortante y portico espacial modelizado co- 
mo un sistema completo (10 grados de libertad) y simplificado ( dos grados de libertad). 



Fig. 2.16 - Modelo dindmico con grados de libertas de rotacion. 
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Ecuaciones de movimiento 

Las ecuaciones de movimiento son las expresiones matematicas que gobiernan la respues- 
ta dinamica de las estructuras. Pueden obtenerse a partir de cualquiera de los principios de 
la mecanica clasica: 

Principio de Hamilton 

*2 '2 

x H = \( E p- E c)dt +\E,dt ; dn H = 0 (2.1) 

h h 

La primera expresion se denomina funcional de Hamilton; donde E p es la energia poten- 
cial, E c es la energia cinetica y E d la correspondiente a fuerzas no conservativas. La segun- 
da expresion permite establecer el equilibrio a traves de una variacion funcional nula. 

Principio de los trabajos virtuales 

Se trabaja en forma similar a lo visto en analisis estatico pero incluyendo las fuerzas de 
inercia y disipativas. 

<Svv, = Sw e (2.2) 


Principio de D’Alembert 

Proporciona el metodo mas directo para obtener las ecuaciones de movimiento de un sis- 
tema dinamico. 

Puede formularse como sigue: “un sistema dinamico esta en equilibrio cuando todas las 
fuerzas que actuan en el mismo, incluidas las de inercia y disipativas, cumplen las ecua¬ 
ciones de equilibrio estatico en cada instante de tiempo. 

Formulacion de la ecuacion de movimiento para un sistema de 1 
GLD 

Tomando el sistema de la figura 2-10, podemos distinguir dos casos: 



Para el modelo (a), aplicando el principio de D’Alembert, tendrfamos: 
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x{t) Fi(t) Fa(t) Fe(t) 

^=if- ◄- ◄- ◄- 

x{t) inercia amort, elastica 


F(t) = Fi(t) + Fa(t) + Fe(t) equilibrio 


Fig. 2.18 - Equilibrio defuerzas para 1 GLD 

A1 aplicar una fuerza exterior F(t), se genera aceleracion, velocidad y desplazamiento para 
un cierto instante “t”; a causa de esto se producen fuerzas: 


i- de inercia 


F,(t) = mx(t) 


(2.3) 


ii- de amortiguamiento 


F a (t) = cx{t) 


(2.4) 


iii- elasticas 


F e (t) = kx(t) 


(2.5) 


equilibrio en el instante “t” 


F(,)-F i (t)-F i (t)-F(t) = 0 
F(t)+F(t)+F(t)=F(t) 
mx+cx+kx=F 


( 2 . 6 ) 

(2.7) 

( 2 . 8 ) 


Se omite (por simplicidad de notacion) la dependencia del tiempo, pero de aqui en adelan- 
te esta se encontrara implicita en toda variable temporal. 


ej: x(t) —> x ; x(t ) —> x ; x(t) —> x 
a(t) —> a ; F(t) —> F ; etc. 
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La ecuacion (2.8) es la de movimiento correspondiente a 1 GLD con carga exterior y 
amortiguamiento. 

Para el modelo (b) de la figura 2-17, el planteo es similar, solo que no tiene fuerza exterior 
aplicada y la fuerza de inercia se ve afectada por la aceleracion total de la masa: 

F i =m[a + x] (2.9) 

entonces, la ecuacion de movimiento queda: 

- m[a +x\-cx — k x = 0 (2.10) 

mx + cx + kx = -ma (2.11) 


La (2.11) es la ecuacion de movimiento para 1 GLD con aceleracion de apoyo (sfsmico) y 
amortiguamiento. Un caso general serfa la inclusion de aceleracion de apoyo y fuerza exte¬ 
rior: 


mx + cx + kx = F-ma 


( 2 . 12 ) 


Formulacion de las ecuaciones de movimiento para modelos con 
multiples GLD 

El modelo con varios grados de libertad mas sencillo es el de edificios de cortante (fig. 2- 
12). Esta basado en que las plantas son infinitamente rfgidas y en que los unicos movi- 
mientos posibles de estas son los desplazamientos horizontales. 

Aplicando el principio de D’Alembert en cada GLD (uno por piso) se obtiene: 


F — F — F — F =0 


(2.13) 


para el piso (r). 

Planteando el equilibrio para todos los GLD, nos queda un sistema de ecuaciones (vecto¬ 
rial) 


F T -F.-F -F =0 

i a e 


(2.14) 
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#Notci: en cidelcinte, los vectores y matrices serdn representados(en general) con minus- 
culas y mayusculas en negrita, respectivamente. 


F ~lI 

F. = -M{x +J a} 
F a - C x 
F e -Kx 


vector de fuerzas extemas 
vector de fuerzas de inercia 
vector de fuerzas disipativas 
vector de fuerzas elasticas 


entonces, el sistema (2.14) puede escribirse: 


M x + C x + Kx=F - M{j a} 

M —» matriz de masas 
C —» matriz de amortiguamiento 
K —> matriz de rigidez 

= [1,1,_11 —> vector con todos sus elementos igual a uno 


(2.15) 


Si bien la (2.15) fue deducida para un modelo de edificio cortante, es una expresion AB¬ 
SOLUT AMENTE GENERAL, inclusive para modelos de elementos finitos, y solo varfan 
las formas de M, C y K. 

Para un modelo sfsmico, la (2.15) se reduce a: 


Mx+Cx + Kx- M {./ a} 


(2.16) 


Y en caso general de porticos 3D, o modelos de elementos finitos, suele sustituirse x por D 
para indicar que cada GLD puede sufrir desplazamientos en 3 direcciones y respectivos 
giros. 


MD + CD + KD-F- M{j a] 


(2.17) 


Nota: en este ultimo caso, J llevara unos en las componentes a las cuales se quiera aplicar 
el acelerograma a, por ejemplo componente x: 


J= [1,0,0,0,0,0,1,0,0,0,0,0,.,1,0,0,0,0,0,.,1,0,0,0,0,0] etc. 
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V y W \_ Capitulo 3 - Respuesta de un Oscilador Simple 

CAPITULO 3 - RESPUESTA DE UN OSCILADOR SIMPLE 


Ecuacion de movimiento y equilibrio dinamico 

Las caracterfsticas dinamicas de un oscilador de 1 GLD pueden estudiarse mediante un 
modelo no amortiguado con vibraciones libres, cuya ecuacion de movimiento es 


mx + kx = 0 


(3.1) 


O X (t ) 


7mmm. 

Fig. 3.1 - Modelo de 1 GLD, no amortiguado. 

La vibracion, del modelo de la fig. 3.1, es inducida por algunas condiciones iniciales, sean 
desplazamiento, velocidad o aceleracion en el instante t = 0. 

Luego, durante las vibraciones no recibe ningun tipo de perturbacion. 

Dividiendo (3.1) por m y usando la notacion: 


2 k k 

CO— — ; co¬ 
ni \l m 

se obtiene: 

x + orx = 0 


(3.2) 

(3.3) 


CO es la “pulsacion” o “frecuencia circular” o simplemente frecuencia del modelo estudia- 
do. Viene expresada en radianes por segundo (1/s). 

La frecuencia ciclica viene dada por 


/ = 


CO 
2 71 


(3.4) 


y se expresa en ciclos por segundo o hertz. 
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Finalmente, otra caracterfstica es el perfodo natural 




La solucion general de la (3.1) o (3.3) puede escribirse: 

x = A sen(cot + V F) 


(3.5) 

(3.6) 

(3.7) 


donde A es la amplitud del movimiento y \\f el angulo de fase. Los valores de A y \\f se cal- 
culan a partir de las condiciones iniciales del problema, por ejemplo para 

*(o )=* 0 ; *(o)=x 0 

resulta: 



Formula de Geiger 

Sustituyendo 



— = X SG : desplazamiento estatico producido 
k 

por el peso G en la direccion del grado de libertad 

entonces: 




(3.9) 
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271 — 

T = X. 


(3.10) 


Utilizando unidades de S.I. la (3.10) queda: 


7 =2.00 X. 


(3.11) 


con X SG expresado en metros para un peso G en Newton. 


Caracteri'sticas dinamicas con amortiguamiento 

El amortiguamiento puede definirse estudiando las vibraciones libres del modelo de la fi¬ 
gure 3-2: 



Fig. 3.2 - Modelo de 1 GLD con amortiguamiento (vibraciones libres) 

Si se toma la ecuacion (2.12) sin cargas ni aceleraciones de apoyo (vibraciones libres) y se 
divide por m se obtiene: 


x + 2 /3x + C0 2 x - 0 


(3.12) 


2/7 = 


(3.13) 


la solucion de (3.12) esta dada en la forma: 


x-e 


(3.14) 


que proporciona la ecuacion caracterfstica: 
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r 2 +2 pr + CO 2 =0 


(3.15) 


Ya que (3.12) es una ecuacion diferencial de segundo orden, lineal, homogenea a coefi- 
cientes constantes. 

Las soluciones de (3.15) son 


r h2 =-j3±^j3 2 -co 2 = 0 


(3.16) 


Segun sea el radicando (3 2 - CO 2 se encuentran tres tipos de amortiguamiento: 


P 2 -C0 2 > 0 -> SUPERCRITICO: la estructura NO VIBRA 


|3 2 -co 2 = 0 —» CRITICO: caso limite (3 = 0 ) —> c r - -2 com 


(3 2 -co 2 < 0 —> SUBCRITICO: la estructura VIBRA con amplitud decreciente 


Este es el caso mas frecuente en ingenierfa civil, por lo que enfatizaremos su 
estudio. Para este caso (subcrftico), la cantidad (|3 2 - CO 2 ) es negativa, lo que 
hace que (3.16) tengaraices complejas: 


r x 2 = ~P ± ico 1 - p 1 = 0 


(3.17) 


con i = 


Llamando frecuencia de vibracion amortiguada a: 


se obtiene: 


co - co 1 -g : 

(3.18) 

6,2 = ~P ± ico x 

(3.19) 

6,2 = ± i 0 ) v 

(3.20) 


En las ecuaciones anteriores aparece la magnitud 

^_c_2Pm 
c 2 com 

r 


(3.21) 
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conocida como fraccion de amortiguamiento critico (en estructuras corrientes 
0.02<^<0.06), tambien: 



co 


o 



c 

2m CO 


Volviendo a la resolucion de la ecuacion de movimiento (3.12), escribimos la solucion ge¬ 
neral en la forma: 


x = c x e nt + c 2 e r ~ (3.22) 

sustituyendo r! y r 2 por la expresion 3.20, se obtiene 

x = Ae~** sen(ffl v f + '¥) (3.23) 

Las constantes A y \\r se obtienen de las condiciones iniciales del problema. 



Fig. 3.3 - Vibraciones libres amortiguadas /1 GLD 


La evaluacion del amortiguamiento en una estructura es un problema esencial en la dina- 
mica estructural. 

El origen de las fuerzas de amortiguamiento se debe a diferentes causas: 

- Rozamiento entre superficies de deslizamiento, en donde la fuerza de amortiguamiento 
es proporcional a la fuerza normal y al coeficiente de rozamiento (hipotesis de Cou¬ 
lomb) 

- Amortiguamiento debido a fuerzas aero o hidrodindmicas 

- Debido africcion interna del material de la estructura 

Generalmente, en el calculo dinamico de estructuras, se utiliza un modelo de gran simpli- 
cidad que caracteriza el amortiguamiento de toda la estructura. Este modelo denominado 
de amortiguamiento viscoso se debe a Kelvin-Voigt y es proporcional a la velocidad. 

F a =-cx 
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Los reglamentos brindan los coeficientes de amortiguamiento para cada tipo de estructura, 
pero puede obtenerse en forma experimental y con un metodo relativamente simple: 


Determinacion practica de 4 

Para amortiguamientos bajos (del orden del 10% de crftico) la relacion entre dos picos su- 
cesivos puede aproximarse: 


x m dM) _ Ae~^sen(co v t n + X P) 
X ,ndx +!) Ae~^ Mn+1 sen{o) v t n+l + *P) 


(3.24) 


pero t n+x =t n + T 


con T =271 CO 


f , = t + 

n +1 n 


entonces 


idx ) 

:(« + l) 


e ^ cotn sen{co v t n + X P) 


-&>\ f n+ 


sen Q) v t n + co v -1- *P 

CO, 


(3.25) 


e ^°* n sen{co v t n +'P) 


•dx ( n + 1 ) 


e ^ a> ‘ n e sen(co v t n + In + 'P) 


X mdxM _ 1 


[n + 1) e 


(O, 


(3.26) 


tomando logaritmo natural: 


In Xmd ^, = ln[/ 2 ™^' 

X mdx\ n + l ) 


(3.27) 


# notar que para amortiguamientos del orden de 
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sA fW~\ - 

£ = 0,1 —> co v = co \-g : 


co -co- 0,995 -» co ~ co 

v 7 v 


Entonces 



x . in) 

max V / 

x . in +1) 

max V / J 


= 


Para el caso de lecturas separadas por N ciclos: 

M/ggfa (* + *)] 
^ 2 7iN 


(3.28) 


(3.29) 


Excitacion periodica 

En la figura 3-4 pueden observarse diversas funciones de carga. De estas, nos interesan por 
ahora, las periodicas y mas particularmente las excitaciones armonicas ya que mediante 
series de Fourier cualquier excitacion periodica puede llevarse a una suma de armonicas 
simples. 



Fig. 3.4 - Tipos de cargas dindmicas. (a) armonica; (b) periodicas; (c) cuasi periodicas; 
(d), (e)fuerzas impulsivas; (f) carga dinamica general; (g) aceleracion sismica del terre- 

no. 
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Excitacion armonica 

Si la carga es de tipo: 


P = P 0 sen (Clt) 


Entonces, la ecuacion de movimiento sera: 


mx + cx + kx = P Q sen(Q.t ) 


(3.30) 


Q.: frecuencia de la excitacion 


la (3.30) puede escribirse tambien: 


x + 2^(0 x + ccrx = — sen(flt) 
m 


(3.31) 


La solucion general de esta ecuacion viene dada por 


X = X, + X 

g h P 


(3.32) 


x h = q e 1 + c 2 e 


(3.22) 


Es la solucion de la ecuacion diferencial homogenea. 


x h =c i e-^e M '‘ +c 2 e-^e- iwj 


(3.33) 


Utilizando matematica para numeros complejos, esta ultima ecuacion puede escribirse: 


x h =e (c x 'sen(O v t + c 2 'cosco v t) 


(3.34) 


x p en la (3.32) es la solucion particular y lleva la forma 


x = A sen fit + BcosClt 


(3.35) 
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Derivando y reemplazando en (3.30) se obtienen las constantes A y B: 


Denotando y = — 

(3.36) 

CO 

P o 1 - y 2 

k (i-fY+(2^rY 

(3.37) 

B _p» -2 4r 

k (i- r 2 ) 2 +(2^rf 

(3.38) 


Condiciones iniciales 

Para el caso en que x(0) = x 0 ; i(0) = i 0 

Pueden calcularse las partes correspondientes a la solucion particular (para t = 0) 


r _10_ _ ^>1 

" m ~k (i-fy + (2trY 

k (i-rf + (24rf 


(3.39) 

(3.40) 


Basados en estas y en (3.32) y (3.34) podemos plantear las siguientes ecuaciones: 

* g (0) = *o = c 2 +x P o (3.41) 

X g (0) = x 0 = -%cg x 0 + %6) x p0 + 0) V C X '+x p0 (3.42) 


y despejarse Ci’ y c 2 ’: 

(3.43) 

(3.44) 
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sA fW~\ - 

Forma de operar : 

1- Se calculan x p0 y x p0 con las 3.39 y 3.40 

2- Se computan cf y c 2 ’ con las 3.43 y 3.44 

3- Se reemplaza todo en 3.32 y se tiene la respuesta en desplazamientos 

4- Para obtener velocidades y aceleraciones se deriva la 3.32 


Nota: la parte de la solucion correspondiente a la ecuacion homogenea incluye el coefi- 
ciente e' sCOt que es una funcion decreciente con el tiempo, por lo que se la denomina transi- 
toria y desaparece antes o despues segun sea el valor de t,. 

La estructura sigue vibrando con una frecuencia practicamente igual a la de la excitacion 
Cl. Esta parte de la solucion se denomina respuesta en regimen. 


Excitacion arbitraria. Integral de Duhamel 


Respuesta a un impulso elemental 


Dada una carga impulsiva como la de la fig. 3-4(d), para condiciones iniciales nulas, pue- 
de calcularse la velocidad y posicion al finalizar el impulso y a partir de ese instante se 
trata de un problema de vibraciones libres. 



Xj = Xj + xdz 


Fig. 3.5 - Respuesta a un impulso 
; x = P(t) m 


entonces: 




m 


(a) 
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Xj = x t + xdT + ^x(drf 


x f =\x(dzf 


(b) 


(b) 


Si dr &s muy pequeno, entonces (b) es despreciable frente a (a), es decir, luego de actuar 
un impulso la masa queda con velocidad pero practicamente no se desplaza. La velocidad 
dada por (a) es la condicion inicial para las vibraciones libres siguientes. 


Reemplazando esta en la (3.23), para un tiempo infinitesimal, se obtiene: 


dx{t) = e-^ (t ~ T) 


dr sen(<y v (t - r)) 

mco v 


(3.45) 


Si el sistema es lineal, pueden superponerse las respuestas de una sucesion de impulsos 
infinitesimales hasta el tiempo generico “t”. La solucion dada por esta superposicion es 
conocida como Integral de Duhamel 


x{t)=^—\p{r)e ^ 0}{t z) sen[co v {t - z)}iT (3.46) 

mco v o 

En el caso de excitacion sfsmica, se transforma en: 

*(?) = —— f a{T)e~^ (t ~ T) sen[d) v {t - t)\It (3.47) 

®v! 

Para el caso en que el problema tenga condiciones iniciales i 0 ^ 0 y x 0 ± 0, la solucion se 
obtiene con el mismo procedimiento dado para las cargas armonicas, solo que la solucion 
particular vendra dada por (3.46) o (3.47). 

Soluciones a la integral de Duhamel : 

(a) Para ciertas funciones de P(x) se encuentran tabuladas las primitivas de (3.46) y (3.47), 
asf como el maximo factor dinamico (FAD) y el instante de maxima respuesta. 
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(b) Para funciones aproximadas por “trozos rectos”, tambien se pueden encontrar solucio- 
nes “exactas” a dicha integral (punto 5.3.4 / pag. 154 / ref. [1]) 

(c) En el caso mas general se puede “integrar numericamente” utilizando alguno de los 
metodos conocidos (Trapecios, Simpson, etc.) 


Factor de amplificacion dinamica 


Es el coeficiente entre el maximo desplazamiento dinamico y el que produciria la carga en 
forma estatica: 


x(t) - 

fad = v /imx 


(3.48) 


Espectros sismicos de respuesta 

Puede definirse el espectro d respuesta (para un acelerograma dado) como los maximos 
valores de la respuesta de un sistema, expresado en funcion del periodo propio. Esto es 
muy util para el diseno de estructuras donde solo los maximos son necesarios. 

Para un sistema de 1 GLD: 

mx + cx + kx - -ma(t) (2.11) 


La solucion a esta ecuacion viene dada por la ec. (3.47) y el maximo valor absoluto sera 

s d . 

Reemplazando S d =|x(t)| . 


Puede obtenerse, para 1 GLD, la expresion exacta del espectro de desplazamientos. Deri- 
vando una y dos veces se obtienen las expresiones para velocidad y aceleracion: 



max 



max 


(Ver ref. [1], cap. 5) 
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Fig. 3.6 - Espectro sismico de respuesta para desplazamientos. Sj 


No obstante pueden hacerse ciertas hipotesis que simplifican las expresiones obtenidas. 
Estas nuevas ecuaciones definen los “Seudo espectros de respuesta”. 

Hipotesis simplificativas: 

i- en estructuras civiles, 0,002 < 4 < 0,2, por lo que puede reemplazarse COy por to 

ii- la funcion coseno que aparece en el espectro de velocidades (al derivar 3.47) puede 

reemplazarse, sin que ello implique grandes variaciones en los valores maximos, por 
una tipo seno 

entonces: 

s d (cd,4)= - — f a(r)e~^~ T) sen[co{t - r)]c/r (3.49) 

S v {(0,<^)= -\a{r)e~^ [t ~ T) sen[(o{t - z)\h (3.50) 

0 max 
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S a i^) = 


co\a{r)e ^ T ^sen[co{t 


-)\h 


(3.51) 


I u I max 

definen los seudo espectros para desplazamientos, velocidades y aceleraciones. 
Es importante notar que: 


(3.52) 

s a =a 2 S d (3.53) 


Por lo que un espectro completo puede representarse en un grafico tipo logaritmico como 
el de la figura 3-7. 



ESPECTRO DE RESPUESTA (AMORTIGUAMIENTO 5 %) 
A— Presa de Pacoima 

B- Presa de Pacoima, roca de cimentaci6n 
C- El Centro 

ESPECTRO DE NEWMARK Y HALL 
D- Eldstico 

E— Elasto-pldstico (/*= 2) 

F- Movimiento del suelo 

DISEfiO SEGON UBC 

G corresponde a R=8 
H corresponde a R=12 


Fig. 3.7 - Comparacion entre los espectros suavizados de Newmarky Hall, los espectros 
de diseno del UBC y los espectros elasticos para los registros de EL Centro y de la presa 

de Pacoima. 


La normativa argentina INPRES-CIRSOC 103 establece la accion sfsmica como una fun- 
cion del seudo espectro de respuesta de aceleraciones S a . Es por eso que solo aparecen gra- 
ficos del tipo de la fig. 3-9, que varfan segun el tipo de suelo y la zona sismica considera- 
da. 
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T(seg) 


Fig. 3.8 - Espectro de diseno suavizado de Seed e Idriss (1982). 


ACELEROGRAMA SINTETICO - 1 





Fig. 3.9 - Seudo espectro y espectro de aceleraciones; (a) corresponde al acelerograma 

(1). (b) corresponde al acelerograma (2). 
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Integration numerica de la ecuacion de movimiento 

Se trata de resolver en forma discreta la ecuacion de movimiento para 1 GLD dada por 
(2.8) y (2.11) utilizando ecuaciones en diferencias que permitan obtener x, x y x en un 
instante t+At en funcion, unicamente, de los valores en el instante t. 

Existen varios metodos para plantear las ecuaciones en diferencias antes mencionadas, pe- 
ro solo nos dedicaremos en profundidad al de Newmark, dado que es uno de los mas di- 
fundidos en software de analisis dinamico mediante FEM. 


Metodo de Newmark (1 GLD) 

4/v- ©-a— | 


i.aft) 

T WJm/M 


mx + cx + kx-F 
con F = f{t) + [- m a(t)\ 


Fig. -3.10 



At = f. +1 -t 

j+1 i 


haciendo el cambio de variable 

T = t-t i ; t < At 

puede suponerse que el valor de 
la aceleracion de respuesta en 


x{r) = Xi +f{r\x M -xi) 

0 para t = 0 
1 para x = At 

integrando (3.55) se obtiene la expresion de la velocidad 


f(T) = 



(3.54) 


(3.55) 


x(t) = X[ + x t T + (x i+l - v.)} f(r)dr 


(3.56) 
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(3.57) 


g(r) = \fWdT 

0 

At 

7 At = j f(t)dt (3.58) 

0 

La expresion (3.56) queda: 

x(t) = Xi + xyc + (x ;+1 - x i )g(r) (3.59) 

que para x=At resulta 

7+i = *i + *i At + feu - 7 )r At 

x i+ i=x i + [(\-y)x i + yx i+l ] A t 

Para calcular los desplazamientos, se integra (3.59) y se obtiene 

T~ T 

x(t) = x t + x t z + X i — + fe +1 - x { )J g(r)dr (3.62) 

Z 0 

particularizando para x=At y llamando 

T 

f3M- =\g(t)dT (3.63) 

0 

se obtiene la relacion final en diferencias, propuesta por Newmark 

71 3 1 , 

7+i = x t +x t T+ --P Xi + px M At (3.64) 

\ z J 

Las ecuaciones (3.61) y (3.64) juntamente con la ecuacion de movimiento (3.54) permiten 
obtener aceleraciones, velocidades y desplazamientos en un instante t+At con solo conocer 
estas magnitudes en el instante t. 


(3.60) 

(3.61) 
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Los parametros y y |3 surgen de un analisis de estabilidad numerica que escapa a este curso 
(ver. ref. [1], [2] y [4]) pero pueden tomarse: 

7 ^ 0,5 
(3 - 0,25 


Desarrollo y forma operativa 

La ecuacion de movimiento (3.54) valuada en t=t i+ i toma la forma: 

mx M + cx M +kx M = F m (3.65) 


Reemplazando en esta los valores en diferencias dados por (3.61) y (3.64), agrupando los 
valores que permanecen constantes en cantidades denominadas por a 0 ,..., a 7 ; se obtiene el 
algoritmo de integracion paso a paso. 


A) calculos iniciales 

1 - determinar las propiedades de la estructura: m, c, k 

2- inicializar x 0 , x 0 y x 0 utilizando, de ser necesario, la ecuacion (3.54) 

3- seleccionar el paso de tiempo y los parametros y y |3 

At; y; (3 At < T/10 ; y=0,5 ; P=0,25 

4- calcular las constantes: 


a o = 


/A r 


a x = 


7 


(3 At 


a 2 = 


1 7 

a= -1 ; a, - -1 ; 

3 2 0 4 p 

a 6 = At(l - 7 ) ; 

5- formar la “rigidez efectiva” 


At 


P At 




a 5 = 


7 

P 


-2 


a 1 =7 At 


k, =k + a 0 m + a^c 


B) para cada paso de integracion 

1- formar el termino de “carga efectiva” en t+At 

f /+ , = F i+ 1 + m(a 0 Xj + a 2 x i + a 3 x i ) + c(a x x i + a A i- + a 5 x t ) 


2- resolver el desplazamiento en t+At 
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= f, 


1+1 


3- calcular aceleraciones y velocidades en t+At 


Capitulo 4 - Respuesta Dinamica de una Estructura con Multiples Grados de Libertad 


sA/W^V 

CAPITULO 4 - RESPUESTA DINAMICA DE UNA ESTRUCTURA 

CON MULTIPLES GRADOS DE LIBERTAD 

Ecuaciones de movimiento y equilibrio dinamico 

Como se vio en el Capitulo 2 un modelo estructural dinamico esta dado por 

- discretizacion espacial (por ej. masas concentradas) 

- discretizacion temporal (por ej. Newmark) 

En el Capitulo 3 se estudio el modelo discreto mas sencillo (1 GLD) mientras que en la 
UT2 se plantearon (sin resolver) las ecuaciones de equilibrio dinamico para sistemas de 
multiples GLD, utilizando el principio de D’Alembert. 

En general, el sistema de ecuaciones diferenciales es del tipo 

MD + CD + KD = F{t ) (4.1) 

En particular, para el caso sismico: 

MD + CD + KD = -M{j a} (4.2) 


Recordar que J es un vector con 1 en la posicion del GLD en el que actua la aceleracion de 
apoyo a(t). 


Vibraciones libres 
Caracteristicas dinamicas 

El sistema que gobiema las vibraciones libres en un sistema de multiples GLD es 


MD+CD+KD=0 


(4.3) 


Como ya se vio, el cambio de frecuencia propia debido a considerar o no el amortigua- 
miento, no es relevante para estructuras civiles corrientes, por lo que seguiremos el anali- 
sis con el sistema simplificado: 
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MD+KD =0 


(4.4) 


Una solucion a este sistema puede ser de la forma 

D =Asen(cot + X F) (4.5) 

A: vector que contiene las amplitudes de las vibraciones 
x f / : angulo de fase 

D(t) = -co 2 A sen(cot+ *¥) (4.6) 

reemplazada en (4.4) queda: 

- M CO 1 A sen(cot + x f / )+ K A sen{cot + x f / ) = 0 

{K-co 2 M}A sen(cot+ x ¥) = 0 (4.7) 

para que haya vibraciones, co^O, por lo que 

sen (cot + x f / )^ 0 

y puede eliminarse de (4.7), quedando finalmente 

{k~0 ) 2 M}A= 0 (4.8) 

Nos interesan las soluciones de A distintas de la trivial. La ecuacion (su resolucion) 4.8 
representa un problema de autovectores y autovalores, en donde: 

K-co 2 M | = 0 (4.9) 

Este determinante puede desarrollarse en la forma polinomica: 

co ln +a x co 2n ~ 2 + a 2 co 2n ~ 4 + ...a n _ x co 2 +a n = 0 (4.10) 

obteniendose la ecuacion caracterfstica. 
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Capitulo 4 - Respuesta Dinamica de una Estructura con Multiples Grados de Libertad 


Cuando K y M son definidas positivas (caso usual en estructuras civiles), de la ecuacion 
caracterfstica se obtienen “n” soluciones positivas cy 2 y en consecuencia “n” valores C0i, 
siendo “n” el numero de GLD de la estructura. 


Los valores cy se denominan “frecuencias propias” o “pulsaciones” de la estructura y los 
“n” periodos propios se calculan 


T i= — i = lv.n (4.11) 

co 

i 

Siendo T : el perfodo correspondiente co 1? que es la frecuencia de menor valor, este se de- 
nomina “periodo fundamental” del sistema. Reemplazando cada 0} en la (4.8) se obtiene 
el correspondiente que se denomina vector de “forma” modal o simplemente “modo”. 
Este vector “modo” contiene la forma que tomara la estructura en cada vibracion. Para ca¬ 
da frecuencia vibratoria, la forma modal sera diferente. Es importante notar que (4.8) se 
cumple para cualquier multiplo de A i; por lo tanto inferimos que no interesa la magnitud 
de las componentes de A- x sino la relacion existente entre ellas. Por esto es conveniente 
“normalizar” los distintos A,. 


Normalizacion de los modos 

1. Por la maxima componente 


dado A = 


a lt 


a. 


a : „ 


A 

9 ,= — 


a , 

imax 


2. Segun la matriz de masas 


puede dividirse 
cada elemento por 

C1 ajmax 


9,1 


(4.12) 


<Pi = 


9,„ 


modo i normalizado 


n=A(M, *)' 


(4.13) 
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Obtencion de los grados de libertad dinamicos 

No siempre todos los grados de libertad estaticos que definen el comportamiento de la es- 
tructura tienen asociada “masa”, es decir que no son necesarios en el analisis dinamico. 
Por ejemplo los porticos de las figuras 2-15 (f) y (h) que son analizados con modelos di¬ 
namicos de menor orden que los respectivos estaticos (GLE). Lo mismo se ve en la figura 
2-16. La forma mas simple de reducir en numero de grados de libertad sin perder precision 
es mediante la: 

Condensation estatica de la matriz de rigidez 

Dada una estructura y su modelo estatico puede subdividirse la matriz de rigidez (asf como 
el vector desplazamientos) de manera de separar las ecuaciones que tienen asociada masa 
(GLD) y las que no la tienen (GLE). 



modelo estatico modelo dinamico 

4x4 2x2 


Fig. 4.1 - Condensation estatica 

Para el ejemplo de la figura 4-1, tenemos que, observando el vector desplazamientos esta¬ 
ticos D s podemos determinar que parte queremos eliminar y cuales grados de libertad de- 
seamos conservar como GLD. 
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reordenando 

{K u -K l 2 K n 'K a )u = F 
KD = F 

donde K es la matriz de rigidez condensada y vale 


(4.21) 

(4.22) 
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(4.23) 


K =K n -K n K 22 K 21 

Finalmente el sistema dinamico de la fig. 4-1 estara representado matematicamente por 


KD + MD = 0 

i 4. 

2x2 2x2 



(4.24) 


Otra forma de obtener la matriz de rigidez del modelo dinamico es invirtiendo la matriz de 
flexibilidad de los GLE que tienen asociada masa. 

Este metodo, si bien menos formal, puede ser de mas sencilla aplicacion para procedi- 
mientos manuales. Para nuestro ejemplo: 



Fig. 4-2 - Obtencidn de K mediante F 1 


Matriz de amortiguamiento 

Dado el sistema 
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MD + CD + KD = F(t) 


(4.1) 


solo nos resta conocer (para encarar la resolution) la matriz C, puesto que en el apartado 
anterior se mostraron dos maneras de obtener K y sabemos (porque usamos un modelo 
discreto de masas concentradas) que: 

0 " 

m, : (4.25) 

L° ••• n \_ 

Conviene definir matrices de amortiguamiento ortogonales, pero para esto debemos, pri- 
mero incursionar un poco en las condiciones de ortogonalidad de los modos propios. 

Condiciones de ortogonalidad 

Las frecuencias y modos propios pueden ordenarse en matrices denominadas espectrales y 
modales, que son, respectivamente: 

co l 0 

n= (4.26) 

_ 0 co N _ 


m 1 

m= : 


Puede demostrarse (verref. [1], [2], [3], [4]) que 



(4.27) 


Esta condicion para los modos puede extenderse a 

(p* M (p j = 0 i^j (4.28) 
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<P, T K<Pj= 0 i*j 


(4.29) 


Si los vectores $ fueron normalizados segun (4.13) las condiciones de ortogonalidad y 
normalizacion pueden expresarse como una unica condicion de ortononnalidad 


® T M® = I 
® T K® = K* 


(4.30) 

(4.31) 


donde I es la matriz identidad y K* es una matriz diagonal 


K*= 


; k t * = (p i K <p. 


(4.32) 


Matrices de amortiguamiento ortogonales 

Una de las hipotesis para lograr una representacion numerica del amortiguamiento de una 
estructura esta dada por suponer que existe un mecanismo de disipacion uniforme de ener- 
gia. En tal caso puede desarrollarse una matriz de amortiguamiento que cumpla las condi¬ 
ciones de ortogonalidad respecto a la matriz modal. 

El “amortiguamiento proporcional” permite definir una matriz que sea proporcional a la de 
las masas, a la de rigidez o a ambas. 


C = a x M + aK 


(4.33) 


La condicion de ortogonalidad: 


<pC(p= 0 i*j 
tfC<p i =2cog i 


(4.34) 

(4.35) 
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Los coeficientes oq y a 2 se calculan a partir de las ecuaciones anteriores (4.35) y (4.33) 


(pi [a,M + a,K\(p i = 2 co i ^ i 

(4.36) 

a x (plM (p t + ajp'K (p. = lore, 

(4.37) 

a x + a 2 tv 1 : = 2 co.%. 

(4.38) 


Recordar que si la normalizacion de los autovectores fue hecha segun (4.13) y (4.14) el 
producto 


<P-K<p=co, 2 (4.39) 

Esto se debe a que la (4.8) puede expresarse para A = q> x y premultiplicarse por (pi 

tf[K-o>?M\p, = 0 
(p i T K<p i -co i 2 <p i T M(p i = 0 

s_ _ / 


(4.40) 

(4.41) 


Determinacion practica de modos y frecuencias 

Este es, por si mismo, uno de los problemas mas complejos de la dinamica estractural. 
Generalmente no es necesaria la resolucion de todas las frecuencias y sus modos, sino que 
solo interesan las primeras “q” que representan las posibilidades ciertas de vibrar pues ne- 
cesitan menos energia de excitacion. 

Cuando el problema es de pocos GLD es posible resolver el determinante de (4.9) o (4.10) 
obteniendose las “n” C0i y luego los “n” (f)„ donde “n” es el numero de GLD. Para proble¬ 
mas con muchos GLD se utilizan tecnicas numericas (aproximadas) para obtener los pri- 
meros “q” pares de (tOi; (f>). En la mayoria de los casos practicos de ingenierfa q « n. 

En las referencias bibliograficas pueden encontrarse descripciones detalladas de varios 
metodos. Nosotros desarrollaremos solo uno que es relativamente simple y preciso. 

Metodo de Stodola-Vianello 

Partiendo de la (4.8) y recordando que se cumple para cualquier multiplo del autovector A 
podemos escribir 
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KA-coMA = 0 

(4.42) 

KA = co MA 

(4.43) 

Suponiendo un vector inicial A 0 conocido (cercano a A) 


K A 1 - co 2 M A 0 

(4.44) 

A 1 ~ K 2 M A 0 

(4.45) 

que lleva a la formula de recurrencia 


A i+1 * K M A 1 

(4.46) 

Nota: una description detallada de este metodo puede encontrarse en f4 /, parte 11.2 
“Vector Iteration Methods”, as1 comoprueba de la convergencia del metodo. 

Procedimiento de calculo 


1- Se propone un vector A 0 inicial. Conviene que los valores correspondan a una defor- 
mada suave, no obstante esto solo acelera la convergencia. 

2- Se calcula 


A 1 = K A M A° 

(4.47) 

3- Se normaliza A 1 


A r = A'{A Ir MA I } _K 
i=l 

(4.48) 

4- Se mejora la solucion 


A i+1 =KM^ 

(4.49) 

5- Se normaliza la nueva solucion 


A !+1 = A ,+l {A i+l,T M A ,+1 } ^ 

(4.50) 


6- Se calcula la frecuencia correspondiente 



7- Se controla la convergencia 


CO -CO 


<TOL —> se vuelve al paso 4 


(4.52) 


<p x =A i+1 ; co { =co i+1 

Se obtuvo el primer autovalor y su correspondiente autovector. 
Siempre: i —> °° ; A' —> cp x ; 


(4.53) 


Obtencion de modos y frecuencias superiores 

El procedimiento es el mismo que para el modo 1, solamente debemos garantizar que el 
vector de “arranque” sea ortogonal a cp\ 

Se supone que un vector inicial 


A ° = ^, 9 , 

j =1 

que para hacerlo ortogonal a cp x le “restamos” la componente q! cp x 


(4.54) 


A*° = A° 2 — q x cp x 


(4.55) 


Para determinar q L , premultiplicamos la (4.54) por \(p x M] y aplicando las propiedades de 
ortogonalidad nos queda 


(p] M]±q j (p i 


(4.56) 


kma: 

q x = Z± T -- 

tfMcp, 


(4.57) 


y en general 


<pMA ] +{ 

Vt = T ~ ,+1 (4.58) 

<P, Mq> t 

y para los vectores de arranque (ortogonales a los previamente calculados) 


A *° = A ; ° — ^ qj cpj 


(4.59) 
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Procedimiento : 

1. SeestimaA ; 0k 

2. Se calcula q; con (4.58) 

3. Se calcula A,-* 0k con (4.59) 

4. Se sigue con el metodo de Stodola propiamente dicho 


Resolucion de las ecuaciones de movimiento en estructuras con 
multiples GLD 

Descomposicion y superposicion modal 

Si bien desarrollaremos el metodo para el caso sfsmico, todo lo expuesto es valido para el 
caso general. 

La ecuacion que gobierna el comportamiento de una estructura de multiples GLD es la 


MD + CD + KD = -M{ja] 


(4.2) 


Las vibraciones libres no amortiguadas se estudian mediante 


MD + KD -0 


(4.4) 


al cual, segun lo expuesto en puntos anteriores, le corresponden “n” pares de frecuencias y 
“modos” que son solucion del sistema de ecuaciones algebraicas 


[k-co 2 m](P = 0 


(4.8) 


Recordemos que los autovectores (p, son ortogonales respecto a las matrices de masa y ri- 
gidez. Es por esto que pueden formar una base completa para el espacio de los desplaza- 
mientos estructurales, es decir, es posible escribir 

D = f J <p,y,(t) (4.60) 
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Donde y^t) es un escalar funcion del tiempo a determinar, llamado “respuesta generaliza- 
da”. 


Sustituyendo (4.60) en (4.2) obtenemos 


.l/\> V ,/!•('Vp ;•(,). A'Vp Vl/! I/./,. (4.61) 

i =1 i=l i=l 

Premultiplicando (4.61) por un (f >, T cualquiera, se cumple que: 


< p]M^ j ( p i =(p T M<p j =M. * 

i= 1 

(4.62) 

1=1 

(4.63) 

Y si, como es habitual, se esta frente a matrices de amortiguamiento proporcionales y or- 
togonales (ver punto 4.4) 

r]c± f V ]c v c* 

1=1 

(4.64) 

Entonces la (4.61) premultiplicada por un $ cualquiera queda como 
GLD 

una ecuacion de 1 

M. * y. (t) + C. * y. (t) + K. * y. (t) = -(p]M J a(t) 

(4.65) 

Recordando que: 


Pj C q>j =C* = Zco^j 

(4.35) 

(p 1 K (p - K * - co 2 

' J ' J J J 

(4.39) 


Y dividiendo ambos miembros por M * = (pjM (p 
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yj W + 2 gjtOj y s (t) + co)y. (t) = 


cp. M J a(t) 

<Pj M <Pj 


(4.66) 


A1 valor 


9jMl_ 

• v;m 9i 


(4.67) 


se lo denomina “ coeficiente de participacion modal ”. 

Finalmente, el sistema (4.2), al ser “proyectado” segun la base formada por la matriz mo¬ 
dal <P= [cp\,... (Pn\, queda en “n” ecuaciones diferenciales de 1 GLD, del tipo: 

y, (t) + 2£.fi>. y. 0 t ) + co]y. (t) = -Q. a(t) (4.68) 


Esta ecuacion puede resolverse con cualquiera de los metodos vistos en la UT3. Una vez 
obtenidos los “n” y,{t), puede obtenerse la solucion estructural mediante (4.60). 

En general, los modos bajos son los que contienen menos energfa elastica de deformacion 
y por ende los que mas contribuyen a la respuesta estructural. 

Usualmente 


i =1 


q < n 


(4.69) 


Procedimiento de calculo 

1- Se determinan los “q” primeros modos y frecuencias 

2- Se resuelven las “q” ecuaciones (4.68) i= 1,... q 

3- Se obtiene la historia en el tiempo segun (4.69) y sus derivadas 

4- Se obtienen las fuerzas elasticas, de inercia y de amortiguamiento para cada instante 


Integration directa de las ecuaciones de movimiento 

Como su nombre lo indica, este metodo no requiere ninguna transformacion previa de las 
ecuaciones de movimiento. Consiste, basicamente en obtener la solucion en una cierta can- 
tidad (discreta) de pasos de tiempo. Es por esto que tambien es llamada integracion “paso 
a paso”. 
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Como ya adelantamos, existen numerosas variantes y algoritmos para la integracion nume- 
rica de las ecuaciones de movimiento, pero nuestros desarrollos estaran dados en el meto- 
do de Newmark . 

En esencia es igual al planteo dado en el Capitulo 3, solo que en vez de tratarse de una 
ecuacion de 1 GLD se trata de un sistema de “n” GLD. 

Procedimiento de calculo 

A) Iniciales 

1. Ensamblar las matrices M, C y K 

2. Inicializar D 0 , D () y D 0 

3. Seleccionar el paso de tiempo y los parametros y y (3 

At <1^/10 ; Y=0,5 ; (3-0,25 

4. Calcular las constantes a 0 ,.. .a 7 

5. Formar la matriz de rigidez efectiva: 


K — K + a 0 M + a^C 


(4.70) 


B) Para cada paso de tiempo 

1. Formar el termino de “carga efectiva” en t+At (* 1) 

r M = F i+ 1 + M(a 0 D. + a 2 D i + afi i ) 

+ C D, + a 4 D j + a 5 D i ) 

2. Resolver el desplazamiento en t+At 

KD M =f M 

* 1: el termino F i+ 1 — F i + {-M J a{t)} se refiere a fuerzas generalizadas 

3. Calcular aceleraciones y velocidades en t+At 

A + , = a o~ D )-a 2 D. - a 3 1) (4.73) 

D M = D - a 6 D, + ci 7 D m (4.74) 

# Notar que para cada paso de tiempo se debe resolver el sistema (4.72) lo que puede 
hacer excesivo el costo de calculo y almacenamiento. 


(4.71) 

(4.72) 


Respuesta maxima utilizando espectros de respuesta 

Las ecuaciones desacopladas (4.68) pueden resolverse utilizando los espectros de respues¬ 
ta: la maxima aceleracion sera 
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(4.75) 


entonces, el maximo desplazamiento es 


y.(t ) = —'-(S ). 

y i v 7 max ,.2 v « / ? 


I max ^ > 2 

co, 


(4.76) 


Podemos, entonces calcular los maximos desplazamientos (en todos los GLD) para el mo- 
do “j”: 




j 2 V a / j 

co j 


(S ). 

D 1 = A.-p- 

j 2 


(4.77) 


(4.78) 


Donde Aj es el vector de coeficientes de participacion modal correspondientes al modo de 
vibracion “j”. Puesto que el maximo en cada grado de libertad no se produce en el mismo 
instante, la respuesta total maxima no es la suma de los maximos de cada modo!! 


D * YZ> 

max max 


(4.79) 


Una forma muy utilizada (y precisa) de evaluar la respuesta maxima (desplazamientos, 
velocidades, aceleraciones y esfuerzos) como combinacion de los maximos modales es: 


R . = X( RJ - f 


(4.80) 


denominada SRSS 
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CAPITULO 5 - ANALISIS DE CONSTRUCCIONES CON 

EFECTOS DINAMICOS DE VIENTO 


T(z): fuerza aerodinamica de 
empuje, paralela a la direccion 
del flujo 

L(z): fuerza aerodinamica de 
deriva, perpendicular a la direc¬ 
cion del flujo 

Generalmente son dominantes las fuerzas de empuje, pero para ciertas estructuras: 

- bajo amortiguamiento 

J- acusada flexibilidad (poca rigidez) 

- construcciones livianas (poca masa) 

las acciones perpendiculares de “deriva” pueden tener efectos significativos que deben ser 
contemplados en el analisis. 

CIRSOC 102-1 —> procedimientos simplificados 



Validez: 

a) construcciones prismaticas o cuasiprismaticas 

b) primer modo dominante y de forma aproximadamente lineal 

c) periodo fundamental Ti > 1 seg 

d) amortiguamiento ^ < 0,01 

Acciones paralelas a la direccion del viento 

(procedimiento basado en el “factor de rafaga”) 



b) 



' / 

/ 


X 


max 




X max {z) = x{z)+x mwc {z) 

x(z): desplazamiento medio 

x max (z): desplazamiento fluctuante debido a la 
turbulencia variable con el tiempo 
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Se define como factor de rafaga 


entonces 


G(z) = l + 


X'max U 

x(z) 


X max (z) = G(z)-x{z) 


Hipotesis basicas para poder reemplazar la accion dinamica del viento turbulento por un 
procedimiento estatico equivalente. 

a) Comportamiento elastico lineal de la estractura. 

b) El modo fundamental de vibracion es una funcion lineal de la altura. 

c) La contribucion de los modos superiores al primero en respuesta se considera 
despreciable, por lo que G(z) = G = cte. 

d) La velocidad media del viento es promediada sobre intervalos de una hora. 

e) La variacion de la velocidad media del viento varfa segun 

\ 0,0706 

Zo,i 

*0.1 y 

,t) 

V 0 : velocidad basica de diseno (m/seg) 
z 0 ,i: parametro que depende de la rugosidad i 
z 0j i: z 0 ,i para rugosidad I (referencia) 

Procedimiento de calculo 
I) Presiones 

La presion dinamica que incluye el efecto de la turbulencia del viento se determina 
mediante 

q\ = G- c z • c 2 ■ q 0 

G: factor de rafaga 

c z : variacion por rugosidad y altura (art. 5.2.4.2, CIRSOC 102) 
c 2 : factor por cambio de tiempo en velocidad media 
(Tabla 3/pag. 20, CIRSOC 102-1) 
q 0 : presion dinamica basica (art. 5.2.3, CIRSOC 102) 
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G = 1 +1,234 K 


B + r ] 



se calcula mediante tablas y abacos que utilizan 
una serie de parametros auxiliares 


En cada nivel se comparara: 


q’z > q z 


NO —> se adopta q z para los esfuerzos 
SI se adopta q’ z para los esfuerzos 


CIRSOC 102 —> q_ = c z ■ c d • q 0 


II) Aceleraciones 

X max (z)= K ■ <7 (z) (paralela al viento) 

K —> fig. 13, CIRSOC 102-1 

<7 (z) valor medio cuadratico de las aceleraciones (pag. 21) 

X 


III) Verificaciones (paralelas al viento) 



a- Dimensionado estructural 
b- Verificacion de confort con graficos y 

tablas en funcion de X max , X max y T 

c- Verificacion de deformaciones admisibles 

— < V < _ 

500 < max < 350 


Acciones perpendiculares a la direccion del viento 


I) Resonancia 


torbellinos de Bernard - Karman 
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Fuerzas de empuie: 


(a la velocidad crftica) 


se admite distribucion uniforme 

T(z) = T. =0,8 C E G-q cr -d 


C E : coeficiente global de empuje CIRSOC 102 
G: factor de rafaga correspondiente a V cr 


q cr —> con V cr // 


Las fuerzas L(z) y T z obtenidas para la velocidad crftica (correspondientes al periodo per¬ 
pendicular y paralelo respectivamente) se suman de la siguiente manera: 

F(z)= :L(zf+T> 

y se debe comparar con las correspondientes obtenidas con q’ z o q z dadas para la velocidad 
de diseno. 
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Criterios de contort en edificios que oscilan 

El confort de los ocupantes de los edificios de gran altura que soportan la accion del viento 
es un tema de primordial importancia en el diseno. 

La respuesta humana al movimiento oscilatorio de las construcciones abarca una extensa 
gama de reacciones, con efectos psicologicos y fisiologicos tales como mareos, ansiedad, 
molestias visuales o temor, llegando hasta nausea aguda. 

Por el contrario, otras persones sienten placer por el movimiento y la experiencia poco 
usual de estar en una construccion que oscila. Adicionalmente, las oscilaciones excesivas 
producen fisuracion de la tabiquerfa, rotura de vidrios de ventanas y cafda de revestimien- 
tos, lo cual incide negativamente en el valor de una propiedad y su rentabilidad. 

Sin embargo, un edificio resultarfa demasiado costoso si se construye o equipa de modo 
que pueda soportar sin movimientos perceptibles una tormenta con vientos huracanados o 
un fuerte sismo. 

En consecuencia, los movimientos son casi inevitables y el problema del diseno consiste 
en mantenerlos dentro de los limites aceptables para no perturbar el confort y el bienestar 
de los usuarios. Por otra parte, el costo del edificio no debe superar los valores normales 
de una aceptable economia. 

Los factores que pueden producir vibraciones en un edificio son numerosos, tales como 
maquinaria en funcionamiento defectuoso, paso de vehiculos pesados por el lugar, impac- 
tos en rampas, vientos fuertes, sismos, etc. los cuales pueden variar durante la vida de ser- 
vicio de la estructura. 

En general la aceleracion es el parametro predominante para determinar aproximadamente 
la naturaleza de la respuesta humana a las vibraciones. Las curvas de la fig. 5.3.1 grafican 
los limites de confort obtenidos del analisis de un gran numero de edificios altos, indican- 
do las maximas aceleraciones aceptables para diferentes frecuencias, dependiendo del uso 
o destino del edificio. 

Los datos se obtuvieron para las aceleraciones pico de las mas fuertes tormentas ocurridas 
durante un periodo de retorno de mas de 5 anos. 
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fig. 5.3.1 


GRAFICA DE CONFORT 


AMPLITUD DE OSCILACION EN FUNCION DEL PERIODO 
PARA DISTINTOS VALORES DE LA ACELERACION 
EN PORCENTAJE DE g 
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La tabla 5.3.1 a continuacion da idea de la magnitud de las percepciones que se obtuvieron 
en experimentos realizados para diferentes niveles de aceleraciones (segun Khan Parmelec 
y segun Chang). 


PERCEPCION HUMANA 

TABLA 5.3.1 

ACELERACION ii (m/seg 2 ) 

No perceptible 

Levemente perceptible 

Perceptible 

Molesta 

(segun Khan y Parmelec) 

ti < 0,004 g 

0,004 g < ii < 0,0075 g 

0,0075 g < ti < 0,02 g 

0,02 g < ti 


(segun Chang) 

No perceptible 

ti < 0,005 g 

Perceptible 

0,005 g < ti < 0,015 g 

Desagradable 

0,015 g < ti < 0,05 g 

Muy desagradable 

0,05 g < ii < 0,15 g 

Intolerable 

0,15 g < ii 
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sA/W^V 

La sensation subjetiva y el comportamiento humano afectado por las diferentes acelera- 
ciones se indica en la tabla 5.3.2 (segun Yamada y Goto). 




TABLA 5.3.2 


NIVEL DE PERCEPCION HUMANA 

RANGO 

ACELERACION 

EFECTO 


(m/seg 2 ) 


1 

0,05 

La gente no percibe el movimiento. 

2 

0,05-0,10 

Las personas sensibles perciben el movimiento. 
Los objetos colgantes pueden moverse algo. 

3 

0,10-0,25 

La mayorfa de las personas perciben el movimien¬ 
to. La oscilacion puede afectar el trabajo de ofici- 
na. La exposicion de larga duracion puede producir 
malestar. 

4 

0,25 - 0,40 

El trabajo de oficina se vuelve dificil o casi impo- 
sible. Aun se puede caminar. 

5 

0,40 - 0,50 

Se percibe fuertemente el movimiento. Hay difi- 
cultad para caminar normalmente. Las personas de 
pie pueden perder el equilibrio. 

6 

0,50 - 0,70 

No se tolera el movimiento y no se puede caminar. 

7-8 

>0,85 

Los objetos caen y pueden lastimar a las personas. 
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GRADUACION 

ACEPT ABILID AD 

ESTADO 


DEL EFECTO 

EN 

PREVISIBLE 


SOBRE LAS 

EDIFICIOS 

DE LA 


PERSONAS 


ESTRUCTURA 

A 


NO 

COLAPSO 

INTOLERABLE 



B 


NO 

DANOS 

MUY POCO 


LOCALES 

TOLERABLE 



C 


SITUACION 

FORMACION 

DEMASIADO 

LIMITE 

DE 

PERCEPTIBLE 


GRIETAS 

D 


EN TAREAS 

FORMACION 

MUY 

INDUSTRIALES 

DE 

PERCEPTIBLE 

PESADAS 

FISURAS 

E 


EN PERIODOS 

SIN INFLUENCIA 

PERCEPTIBLE 

BREVES EN 

EN ESTRUCTURAS 



VIVIENDAS 

CORRIENTES 

F 


PERIODOS 

SIN 

ESCASAMENTE 

LARGOS 

INFLUENCIAS 

PERCEPTIBLE 

EN VIVIENDAS 


G 


PERIODOS 

SIN 

NO 

LARGOS 

INFLUENCIAS 

PERCEPTIBLE 

EN VIVIENDAS 



CUADRO AD JUNTO A GRAFICA DE LA fig. 5.3.2 
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Capitulo 6 - Analisis de Construcciones con Efectos Sismicos 


CAPITULO 6 - ANALISIS DE CONSTRUCCIONES CON 

EFECTOS SISMICOS 


Definicion numerica de la accion sismica 

La definicion correcta de la accion sismica es un problema al que se le debe dar la mayor 
importancia en un analisis estructural sismico. 

La solucion a este problema parte, generalmente, de los datos experimentales proporcio- 
nados por la sismologia. 


Parametros de Prediccion de las 

terremotos pasados -> caracteristicas sismicas 

en una region de los terremotos que 

afectaran dicha region 


Tradicionalmente, la fuerza destructiva de un terremoto ha sido expresada en funcion de la 
aceleracion maxima del terreno, pero se han observado danos moderados a aceleraciones 
muy altas. 

Por lo tanto, se deben tener en cuenta otras caracteristicas: 

- intensidad, duracion 

- contenidos de frecuencias 

- secuencia de choques, etc. 

Por ejemplo, se han utilizado, ademas de los valores maximos de la aceleracion, velocidad 
y desplazamiento del terreno, el espectro de amplitudes de Fourier, el espectro de seudove- 
locidades, intensidad espectral (Housner), el valor medio cuadratico de las aceleraciones 
en la fase “fuerte” de un acelerograma, etc. 

Ademas, la accion definida estara totalmente relacionada al tipo de analisis estructural que 
se va a realizar: 


Analisis lineal 

desacoplamiento 
-» -> 
modal 

espectros de 

respuesta 




Analisis no lineal 

integracion 

-» -» 
directa 

acelerogramas 
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# Tambien existen definiciones de acciones sfsmicas a partir de la teorfa de procesos esto- 
casticos, en donde la respuesta estractural en si se obtiene probabilfsticamente. 


No profundizaremos sobre estos temas, pues escapan a nuestro curso. (ver ref. [1 y 2]) 


Definicion mediante espectros de respuesta 

Es la forma mas usual de definir una accion sfsmica, dado que se obtienen descripciones 
de las caracterfsticas mas importantes de la respuesta estructural sin la necesidad de dispo- 
ner de una historia en el tiempo de la excitacion y la respuesta. 

Otra ventaja es que un espectro puede “modificarse” en base a las caracterfsticas del lugar 
de emplazamiento de la estructura, sin necesidad de conocer los detalles de la excitacion. 


Espectro 
de respuesta 




representaciones graficas de valores 
aproximados de la respuesta maxima 
de un sistema de 1 GLD lineal y elastico: 


x+ 2^co ■ x+ or • x = - a(t ) 

S , (a, co, £) = - — • f a(r) ■ e ~^ o (t ~ T) - sen co (t -t)- dr 
co J 0 

S v (a,6J,g)=6J- S d (a,6J,g) 
S a (a,co^)=co 2 ■S d {a,o ,£) 


( 2 . 11 ) 

(3.49) 

(3.52) 

(3.53) 


Puede realizarse un analisis de las variaciones de los valores de las curvas S d , S v , S a , en 
funcion de las caracterfsticas de la estructura. 


Recordar que CO = 

m 

Para frecuencias propias altas (en comparacion con la del movimiento del terreno), la 
(2.11) quedarfa (despreciando los dos primeros terminos): 


estructuras 
muy rfgidas 


CO 2 ■ x ~ -a(t ) 

< .\si<y^°o;S a ^a max 
(copia la aceleracion maxima del terreno) 
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Para frecuencias bajas, por el contrario, la (2.11) queda: 


estructuras 
muy flexibles 


x ~ -a(t) 

si<y^0;S d d max 

(copia los desplazamientos maximos del terreno) 


Para frecuencias intermedias, se produce una amplificacion del movimiento del terreno en 
su paso a traves del filtro estructural. 


PM'T)=—-, PMJ)=—, /3A^t)= 4 a - 


(6.1) 


Las condiciones locales del terreno, tales como grosor y propiedades de los estratos que se 
encuentran entre la roca firme y la estructura modifican los espectros de respuesta. 



Parte 1 / Cap. 7 
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(T 7 } 

V 

73 

b • 

2 



[Yj 



T[ <T<T 2 


t 2 <t 


a s ,b,T x y T 2 —> tabla 4 —> / 


^zona sfsmica^ 
v suelo y 


Para casos con 0,005 < 4 < 0,05 (lrnea a trazos) 



T < Tj 


Tj < T < T 2 


t 2 <t 


/ A : factor de amplificacion por amortiguamiento 


f A 


5 

t, (expresado en %) 


Definicion mediante acelerogramas 

a) Acelerogramas reales 

Basar un calculo sismico en uno o varios registros disponibles en una zona implica 
un alto riesgo de definicion incorrecta de la accion. 

b) Acelerogramas sinteticos 

Tiene grandes ventajas respecto al anterior: 

posibilidad de generar senales de corta duracion que tengan el mismo efecto en 
las estracturas que el del terremoto real que se quiere simular => economfa 
computacional 

al estar generados en funcion de un espectro de diseno, se tiene en cuenta 
(aprox.) las condiciones locales del suelo 

La generacion de acelerogramas sinteticos requiere procedimientos matematicos basados 
en procesos estocasticos y expansiones en series sinusoidales. 
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Existen muchos procedimientos y algoritmos que permiten generar acelerogramas sinteti- 
cos a partir de un espectro dado. 

INPRES - CIRSOC 103 (Parte 1 - Cap. 14.3.1) 

Los acelerogramas deben satisfacer: 

a) 

a (t)ma, Z a s ■ y d 

a s : ordenada al origen del espectro correspon- 
diente 

y d : factor de riesgo, art. 5.2 

b) El espectro elastico (para t, = 0,05) debera cumplir: 


Area c >y d ■ Area D 

S c a > 0J • • y d 

para todos los puntos 

Cantidad de acelerogramas a emplear: 

Grupo A y B n° > 3 acelerogramas 

Grupo A 0 n° > 4 acelerogramas 

Para diseno y verificaciones se promediaran las envolventes de solicitaciones y deforma- 
ciones obtenidas para cada acelerograma, pero en dicho promedio no se incluiran valores 
inferiores al 85% del maximo encontrado. 
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Metodos de analisis segun INPRES - CIRSOC 103 


Procedimientos con fuerzas estaticas equivalentes 

a) Metodo estatico (art. 14.1) 

b) Procedimiento aproximado (cap. 16) 

Metodos dinamicos 

a) Analisis modal espectral 

b) Superposicion modal paso a paso 

c) Integracion directa paso a paso 

6 ... Metodos dinamicos INPRES - CIRSOC 103 

a) Analisis modal espectral 

• excitacion sfsmica traslacional 

• materiales lineales y elasticos (segun art. 12.1 / pag. 47) 

• (s a ).=s^-r d -g 

S a : del espectro correspondiente 

y d : factor de riesgo 
g: aceleracion de la gravedad 

• se podra considerar capacidad de disipacion de energia por deformaciones 
anelasticas de la estructura 

g R 

las deformaciones calculadas con este criterio deberan ser amplificadas 
multiplicandolas por la ductilidad global p 

para el caso anterior, es imprescindible que la estructura se comporte en 
forma uniforme, de manera de garantizar la ausencia de concentracion de 
rotulas plasticas 

Se debe verificar que el corte (con analisis modal espectral) en cada direccion no resulte 
inferior al 75% del corte obtenido con el metodo estatico; de no verificarse se incrementa- 
ran los efectos del analisis modal de la siguiente manera: 
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R* 


modal 

espectral 


OJ5-V 0est 

modal 

espectral 


■R 


modal 

espectral 


Modos a considerar : q > 3, pero ademas se deberan incluir todos los modos cuya contribu- 
cion al total sea mayor que 5% de la contribucion del modo fundamental. 


Modelo de edificio de cortante 


h 


m. 


w, 


in 


m, 


k i 


x. 


- X 

X 


- X 


r-1 


a 


(0 


(p T ■ M -{l} 

Q,= T „ = v (6-1) 

L m , ■ <Pri 




(=i 

0 " 


<Pu 

con M = 


m r 


; <p,=' 

<Pri 


_0 


m r 


Vni, 


# observar que si en (6.1) se susti- 
tuye la masa m r por el peso w r el 
valor de Qj no cambia 


si tpf-M-p. =1 
Q i =Y J m r -(p ri =Y J w r -(pr 


i—l 


i—l 


Peso efectivo modal 



—> si —> 

r =1 

q => apto 

Desplazamientos modales maximos 


iw, >o,9- ±w r 

i =1 


i =1 


peso total de la 
estructura 


( X ,L, = 9, '-'I = f, Q, ■ 


(sA 


co: 


(6.4) 


(6.5) 
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Desplazamientos modales relativos entre pisos 


A n = Ur ), - Ur-1 ), 


( 6 . 6 ) 


Fuerzas sfsmicas equivalentes (modales) 


F ri = Xri = m r 


(6.7) 


pero 


Xri) =<Pri-Xi(t)„ 


(6.8) 


Como 


x <( t ) max =Q i -( S a) i 

Xn 1 - Vri ' Qi ' CO; 


F- = <Pn ■ Qi ■ {S a l ■ m r = (p t ■ Q t 


fuerza sfsmica equivalente 
en el nivel r para oscilaciones 
en el modo i 


(6.9) 


( 6 . 10 ) 


en funcion del peso w, 
y la fraccion de g 


Vectorialmente: 


M (S ' 

<pj - M- <p : 


( 6 . 11 ) 


¥f=M-g>,-Q r {S,l 


( 6 . 11 ) 
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Cortante modal 


En nivel “r”: 


V. = Y J F S 

n Lu ]i 


( 6 . 12 ) 


En la base: 


o tambien: 


v,=2X s 

r =1 

v =w .^ 

i i 

8 


(6.13a) 


(6.13b) 


Comparacion con INPRES - C1RSOC 103 


V = — . s • w 

m wj ^ am 1 m 

K 


(Formula reglamento) 


y d : factor de riesgo 

R: coeficiente de reduccion por disipacion 


i = m 


(S ) 

V,=f r — - W, 
8 


(Formula 6.13b, multiplicada por un factor/i) 


Distribucion en altura de V m segun INPRES - CIRSOC_103: 


T? _ W k ' |fifon _ y 

km n y m 




m: modo 
k: nivel 

F km : fuerza sfsmica equivalente para el modo m en el nivel k 
Segun la formula (6.10) y (6.13a): 

T/ ^ Qi'( S a)i sr, n 

Vi =L ( Pn ■ Qi ■ W r -=- LVri ’ K 

>'= 1 g g r= 1 

si queremos emplear la notacion 

F s =£.-V. => F s =£ . • - Y<p • w 

g r =1 

comparandola con (6.10) se deduce que 
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Zri = 


(p ri -w r _ W k -(f) km 

i<Pn-W r iw r <t> ir 


notacion INPRES - CIRSOC103 


Momento de vuelco modal (para el nivel r) 


K, = t F iXhj -ft,) 


si r = 0 —> Mn; v = momento de vuelco en la fundacion 


Momento de torsion modal 



M' ri = e r ■ V rj corresp. al piso “r” 


Esfuerzos seccionales modales 

Para cada subestructura vertical con la parte proporcional a su rigidez que le corresponde 
de la fuerza sfsmica equivalente, se calculan los esfuerzos seccionales para cada modo en 
estudio. 


a ■ F 

p n 



subestructura "P" 


nivel "r" 


CC X + OC 2 + ...dp + ... OL m — 1 

m: n° de subestructuras verticales 


(se deben adicionar los efectos torsionales si los hubiese) 


CALCULO DE LA RESPUESTA TOTAL 
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